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Abstract. We use the Star Track stellar population synthesis code to analyze properties of double compact object 
binaries as sources of gravitational waves. Since the distribution of lifetimes of these objects extends up to the 
Hubble time we conclude that a proper calculation of the expected rate must include a full cosmological model. 
We present such model, calculate the expected coalescence rates, and analyze the intrinsic sensitivity of these rates 
to the model assumptions. We find that the rate alone is a very poor estimator of the underlying stellar evolution 
model. However we show that the distribution of observed chirp masses is very sensitive to the underlying stellar 
evolution model, while it is very insensitive to the underlying cosmology, star formation rate history and variation 
of detector sensitivity. 
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1. Introduction 

The large interferometric gravitational wave detectors are 
entering the realm of astronomy: LIGO has now completed 
it first two science runs, and is improving its sensitivity, 
VIRGO preparations are proceeding very well. GEO 600 
has achieved an astonishing duty cycle above 97%, and 
TAMA 300 has been taking data in coincidence with 
LIGO. The sensitivity of all the instruments is improv- 
ing and this leads to increased hopes for detecting grav- 
itational waves. While a detection of gravitational wave 
will be a huge triumph by itself, the following question 
appear: can there be any useful astrophysics done with 
gravitational waves? can gravitational wave observations 
provide any constraints on astrophysical models, and be 
used to answer questions related to what is known as stan- 
dard astronomy? 

Coalescing compact object binaries are one of the most 
promising, if not the most promising, sources of gravi- 
tational waves. First, we know from electromagnetic ob- 
servations that they exist and that they emit gravita- 
tional waves. Moreover, the known sources will coalesce. 
The research on the properties of coalescing compact ob- 
ject binaries has concentrated so far on calculation of 
the expected rates. This problem has been approached 
in two ways. The first approach was based on analyzing 
the properties of the existing systems, and considering all 
the possible selection effects that may affect detectabil- 



ity o f such systems l|Naravan et alJll99ll iKaloeera et alJ 
l200lh . The drawback of this approach is the small num- 
ber statistics, or even zero object statistics, in the case 
of black hole neutron star or double black hole bina- 
ries. Moreover, the estimates of different selection ef- 
fects carry additional uncertainty. A second approach is 
based on studying the stellar evolution processes and de- 
tailed analysis o f the formation paths of double compact 



object binaries lIliEunov. et all Il997bt iFrver et al.l 



Bethe fc Brownl 



1998; 



1998; 



I Portegies Zwart fc Yungelsor 
Fryer et al.' '1999; 'Bclczynski & Bulik '1999'; 'Bufik et alJ 
.1999; ..Bclczynski et al. .2002c .Nutzman et ah, 200^ 
Within this approach the main problem is insufficient 
knowledge of certain processes in the stellar evolution. 
Such processes were parameterized and the results appear 
to vary significantly with some of this parameters. The 
uncertainty on the expected rate of compact object coa- 
lescences is between two and three orders of magnitude. 

In our previous paper jBulik fc BelczvhskilEiol we 

pointed out another possibility: we have shown that mea- 
surement of chirp masses of coalescing systems carries a 
lot of information about the underlying stellar evolution. 
This calculation used a set of simplifying assumptions: 
Euclidean space geometry and a constant star formation 
rate. In this paper we relax these assumptions and con- 
sider a more general model. In section 2 we describe the 
population synthesis model and argue for the need to con- 
sider a fully cosmological model of the source distribu- 
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tion to calculate the distributions of observable quantities. 
Section 3 contains an outline of such calculations, and a 
discussion of the dependence of the final results on vari- 
ous cosmological parameters and assumptions. Finally we 
finish we conclusions in Section 4 

2. Population synthesis implications 

In th is work we use the Star Track population synthesis 
code ( Belczvnski et al.l2nn2cl) . The single star evolution in 
the code is b ased on the a p proxi mate analytical formulas 
compiled by iHurlev et all |20o3). The single star evolu- 
tion description includes such stages as the main sequence, 
Hertzsprung gap evolution, red giant branch, core helium 
burning, asymptotic giant branch, and the helium star 
evolution. The end product of the stellar evolution can be 
a white dwarf, a neutron star, or a black hole. The binary 
evolution description includes such processes as variation 
of the orbits due to wind mass loss, tides, as well as mag- 
netic breaking, and various mass transfer modes: conser- 
vative, quasi dynamic, common envelope (CE) evolution. 
We also take into account the variation of the structure 
of a star in response to accretion (e.g. rejuvenation) and 
a possibility of hypercritical accretion onto compact ob- 
jects. Supernovae explosions are treated in detail usin g 
the results of hydrodynamical simulations llFrveT<ll999l) . 
We take into account both direct and fall back black hole 
formation. Finally, we also investigate various kick veloc- 
ity distributions. The initial mass of the primary is d rawn 
from a power law distribution oc M~^'^ (Sc^d^^^Q), the 
initial mass ratio distribution is flat (KuiDenilQSSj), the 
eccentricities are drawn from a d istribution oc e l)Heggij 
Il975t iDuauennov fc MavMlll99l|) , and the orbital sepa- 
ration distribution is flat in log a llAbtlll983l) . 

Apart fr om the standard model ( model A) described 
in detail by iBelczvnski et al.l ((2002cl) we also consider a 
variety of different population synthesis models where we 
vary parameterizations of stellar evolutionary stages. A 
list of models used later on is presented in Table 1. 

In considering the standard model of population syn- 
thesis (model A) we evolved Ntot = 10^ initial binaries. 
All other models were were run with Ntot = 4 x 10^ ini- 
tial binaries. The number of coalescing compact object 
binaries obtained in each simulation is listed in Table 1. 
An example of an evolutionary path leading to formation 
of a double black hole binary is presented in Figure ^ 
In this particular example, the evolution starts with two 
massive stars (primary mass Mi = 53.0Mq, secondary 
mass Ml = 28.4Mq) at Zero Age Main Sequence on a 
rather wide (semi major axis of orbit a — 176Rq) and ec- 
centric orbit (eccentricity e = 0.4). More massive primary 
evolves first off Main Sequence (MS) and starts evolution 
on a giant branch. Due to a rapid radial expansion, it 
fills its Roche lobe circularizing the orbit and initiates the 
first mass transfer (MT) episode (stage II). This transfer 
is stable: rather small orbital contraction and we assume 
non-conservative evolution: primary looses its entire en- 
velope, half of which is lost from the system and only 



M,[Mj,] STAR 1 STAR 2 MjMj, 



(I) 53.0 



(II) 43.1 



(III) 15.5 



(IV) 8.6 



(V) 6.2 



(VI) 6.2 



(VII) 8.0 



(VIII) 8.0 



(IX) B.O 



a.[Rj e 
176 0.4 
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Fig. 1. An example evolutionary scenario leading to for- 
mation of a double black hole binary. For details see the 
text. 



the remaining half is attached to the companion. After 
MT episode we are left with a bare Helium core of the 
primary and now the more massive rejuvenated MS sec- 
ondary (stage III). Primary, evolves fast and ends its life 
in a SN lb explosion, inducing a slight eccentricity on the 
binary and widening the orbit (IV). As a result of the 
core coUapse/SN event a first BH is formed (V). Now, 
the secondary follows the early path of its companion. 
It leaves MS and becomes a giant star. Radial expansion 
leads to the second MT episode, which due to the extreme 
mass ratio of two components (stage VI) , is dynamically 
unstable. The system goes through a Common Envelope 
(CE) phase, with drastic shrinkage of the orbit, while al- 
most entire secondary envelope is ejected from the system. 
Small fraction of the envelope is accreted onto the BH and 
the secondary becomes a massive helium star (stage VII) . 
Further nuclear burning of elements in the interior of the 
helium star brings it finally to a core collapse and SN lb 
explosion (stage VIII). The second BH is formed, and we 
end up with two massive stellar BHs orbiting each other 
on a tight and eccentric orbit (stage IX). 

For each compact object binary we note the masses 
of individual objects mi and TO2,and the lifetimes tufe 
of the system from its creation on the zero age main se- 
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Fig. 2. Distribution of the lifetimes of the double neu- 
tron star binaries (solid line) , black hole neutron star bina- 
ries (dotted line), and double black hole binaries (dashed 
lines) . 



quence to formation of the q compact object binary and 
then due to the gravitational wave energy loss. Within our 
code the distinction between a neutron star and a black 
hole is based solely on its mass. In our calculations we 
have been assuming a maximum mass of a neutron star 
^'^max — 3 -^^0) all compact objects heavier than M^^^ are 
considered to be black holes. Given that we can distinguish 
three different types of binaries in our calculations: dou- 
ble neutron star binaries (NSNS), black hole neutron star 
binaries (BHNS) and double black hole binaries (BHBH). 
It appears that some properties of these binaries are dif- 
ferent. 

We plot in Figure |3 the distribution of the lifetimes 
of these three types of binaries obtained within model 
A. A majority of double neutron star binaries are rather 
short liv ed, with the lifetimes in the range from 10 to 
SOMyrs l)Belczvnski et all 120024 - The mixed BHNS bi- 
naries lifetimes span w very wide range from roughly 
lOMyrs until the Hubble time with no really preferred 
interval. The BHBH binaries live much longer their life- 
times extend from about lOOMyrs to the Hubble time. 
Thus the currently merging NSNS binaries originate in 
stars formed in relatively recent starbursts, while the 
BHBH and BHNS binaries origin ate in stars that have 
been formed a few billion years agoi'Lip unov et all ( 1997a^ 
and recentlv iBulik fc Belczvhski (20Q2j) have shown that 
the observed sample of coalescing binaries is dominated 
by the BHBH binaries. Since they are so long lived and 
the star formation rate was probably much higher in the 
early Universe the accurate calculation of their observed 
properties requires taking into account a full cosmological 
of distribution and evolution of such sources. 



3. Expected numbers of observable mergers 

The signal to noise from a coalescence of a binary in grav- 
itational wav e interferometers has been calculated by sev- 
eral authors jCher noff fc Finn 19931 iBonazzola fc Marckl 
Il994t [Flanagan fc Hughea.1998^1 . In the frequency region 
where the detectors are most sensitive it is 
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(1) 



where M = (mim2) (mi -I- m2)~ is the chirp mass, 
z is the redshift, d]^ is the luminosity distance, and the Ai 
is the distance at which a coalescence of a 1 Mgdot binary 
produces (s/ N) = 1 in the detector . In a m ore detailed cal- 
culation - see lFlanagan fc HughesI ((l99^ for the cases of 
LIGO I, LIGO II, and LISA - equation Q becomes more 
complex due to the finite size of the frequency interval 
available for a given detector. This formula is valid with 
the accuracy better than 10 percent for chirp masses be- 
low < 20 Mq. In the following calculations we assume a 
flat cosmological model with the matter density fim = 0.3, 
and cosmological constant — 0.7, the Hubble constant 
Hf) is 65km~^Mpc~^. We assume a model of star forma- 
tion history rate presented as a thick solid line in Figure|21 
(^owan- Robinson 1999.) , but we also consider a model in 
which the s tar formation rate ceases at z > 2 - thin line 
in Figure Ol l|Madau et al.lll996(l . and a flat star formation 
rate. 

Finding the sampling redshift zm requires solving 
equation (^), with a required value of (S/N), where we 



insert di, = (1 



z)dprop, and dprop is the proper distance 



,(z) = / {l + z)c 
Jo 



dz . 



Here we have introduced the cosmic time t: 



1 



Ho{l + z)[{n^{l + zf + (1 - r!„)]i/2 



(2) 



(3) 



in the flat space time. 

Formation of the compact object binaries is directly 
connected with the star formation rate. After a binary is 
formed it evolves for a time tuf^ until its merger. This 
time includes the stellar evolutionary time needed to form 
a compact object binary and then the life time of such 
binary due to gravitational wave emission. The coales- 
cence is delayed by tufe with respect to the star forma- 
tion. Let us denote by F{J^,t,z) the formation rate of 
binaries with the chirp mass M and a lifetime i as a 
function of redshift. The star formation rate history in 
the Universe SFR{z) can now be combined with the re- 
sults of the population synthesis code to flnd F{Ai, t, z). 
To this end we need two quantities: the average stellar 
mass (Af»), so that we find the star number formation 
rate S F R{z) / {M^,) , and the fraction of the stellar popu- 
lation that we simulate with the stellar population syn- 
thesis code fsim, to obtain the number formation rate of 
our binaries as a function of z: SFR{z)fsim/{Mf). For 
the assumed slope of the initial mass function a = —2.7, 
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Table 1. Description of different population synthesis models used here. We list the number of binaries produced in 
each simulation. 



Model 


Description 






N produced 


A 


standard model described inlBelczvnski et alJ (|2002(J). 


33437 




hilt witVi Tti MI — 1 ^rivrs 








Bl 








47690 


B7 


ciino'lp A/Tj^ Ywrpl 1 i n Ti wjitVi rr — ^0 Q ^ 
ollltilVj IVJ-dA WCllitXli WlLil c — r^iii o 






47715 


Bll 


single Maxwellian with tj = 500 km s""'^ 






5246 


B13 


IPaczvnskl Jl99[J) kicks with 14 = 600km s~ 


1 




5247 


C 


no Vivnpr — prif ipftl apprptinn nniri IVSi /Rl~r in 


CEs 




8416 


El 


rVf-^-r^ \ — 1 






1003 


E2 


rvriT-i V X — n 






5538 


E3 








17595 


Fl 


iiicLSS fr&ction cLCcrstsdi = 0.1 






4461 


F2 


mfiQQ frftptinn f\ccT(^\(^c\' f — 1 

lllcLOO 11 ClL^LlWll CJbKj\^L C L v_i . Ig^ _1 






9885 


Gl 


wind clicLnged by /wind = 0.5 






12345 


G2 


wind chaneed bv f,,r;nH = 2 






11482 


J 


primary mass: oc Mj~^'^^ 






13903 


LI 


angular momentum of material lost in MT; 


j = 


0.5 


13529 


L2 


angular momentum of material lost in MT: 


j = 


2.0 


3335 


Ml 


initial mass ratio distribution: ^{q) oc 






1363 


M2 


initial mass ratio distribution: oc (f' 






22097 


O 


partial fall back for 5.0 < Mco < 14.0 Mq 






9193 


s 


all systems formed in circular orbits 






8306 


Zl 


metallicity: Z = 0.01 






8520 


Z2 


metallicity: Z = 0.0001 






10481 



the minimal stellar mass of 0.2 AIq, maximal stellar mass 
of 100 M0, and binary fraction of 0.5, the average stellar 
mass is (M*) = 0.87 Mq, and the fraction of stars that we 
simulate is fsim = 1-24 x 10"'^. The numerical estimate of 
the formation rate of binaries in a given interval between 
M and M + dA4 , with lifetimes between t and t + is 



F{M,t,z)dMdt = SFR{z 



(4) 



where Ni is the number of binaries in a simulation with 
the chirp mass and the lifetime in a given interval and 
Ntot is the total number of simulated binaries. 

Binaries coalescing at a given redshift zq originate from 
binaries formed at different earlier times. The rate of co- 
alescences of binaries with a given chirp mass A4 is then 
given by 



dfc. 



dM 



izo,M) = / dt'FiM,t',zf), 



(5) 



where the source formation redshift Zf is obtained by solv- 
ing: 



t' 





dt 


J Zo 


dz 



dz . 



(6) 



We can now proceed to calculation of the observed rate 
of coalescences. This calculation is very simil ar to the ones 
performed in the case o f gamma-ray bursts l|Totanilll999l 
iBelczvnski et al.ll2002 b^. An instrument can detect signal 
characterized by the chirp mass M. if it lies closer than the 



sampling redshift zm The rate at which such instrument 
will detect binary coalescences is 



dR 
dM 



df, 



coal 



dM 



1 dV 

1 + z dz 



dz 



where M' = M{1 + z) ^ and 



dV , d^r 
— = 47rc--— ^ 

dz (1 + z) 



(7) 



(8) 



is the comoving volume element. 

The results of the differential rate calculation are 
shown in Figure0]where we present the model requir- 
ing a signal to n oise (S/N) > 8 and using the sensitivi- 
ties calculated by Flanagan & Hughes fl 998j) . LIGO I and 
LIGO II and bv HeUo 1,1998) for VIRGO, i.e.: Aligo i = 
100 Mpc, AviRGO = 150 Mpc, and Aligo i = 2200 Mpc. 
There two main effects influencing the presented in 
Figure ^ the first is that with increasing sensitivity a de- 
tector detects coalescences from a larger volume and the 
second is due to the fact that the observed quantity is 
the redshifted chirp mass (1 -I- z)M. A highly sensitive 
detector like LIGO 2 shall detect coalescing binaries at 
non negligible rcdshifts and this leads to the spread of the 
distribution of observed chirp masses to higher values. 

At this point we must also note the uncertainties in 
the estimates shown in Figure 0] especially coming from 
the product SFR{z)fsim{M^)~^. The estimates of the star 
formation rates are uncertain by a factor of at least two 
l|Madau et alJll99^ . Both fsim and (A/*) depend strongly 
on the assumed slope of the initial mass function, and the 
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Fig. 3. The star formation r ate history models used here: 
thin l ine l|Madau et alJlQQ^) , thick line l(R,owan-R,obinsonl 
Il999l) 
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Fig. 4. The differential detection rates expected for the 
existing long base interferometers. 

lower mass cutoff of the stellar population, and weakly on 
the binary fraction. The low mass stars, not taken into ac- 
count in our simulations since they do not lead to compact 
object binaries, dominate the stellar population in num- 
bers and in mass. We have assumed a power law shape of 
the initial mass function, yet the absolute numbers of the 
high mass progenitors of compact object binaries will be 
affected by deviation of the initial mass function from a 
simple power law. A deviation of the slope by 1/2 changes 
the number of high mass stars in relation to the total pop- 
ulation by a factor of ten! Let us summarize: In the cal- 



culation of the rate using equation iQ, there are at least 
the following uncertainties: the unknown binary fraction, 
could be as high as a factor of two, the unknown shape of 
the initial mass function leads to uncertainty by a factor of 
up to ten. The differences between the galactic coalescence 
rates obtaine d in the framework of va rious models of stel- 
lar evolution ijBelczvnski et al.1l2002c|) amount to a factor 
of about 30. The difference in the expected rate within 
various stellar evolution models is comparable to the in- 
trinsic uncertainty in the rate. Thus given a measurement 
of a rate, answering an inverse problem, i.e. what does the 
rate tell us about the underlying stellar evolution, shall be 
next to impossible. Here we have assumed that the mea- 
sured rate has no uncertainty, yet we know already that 
the interferometric detectors exhibit a very non stationary 
noise which hampers an honest estimate of the time space 
volume surveyed. 

4. Distribution of chirp masses 

We note that the shape of inferred distributions of chirp 
masses in Figure 0] is nearly identical for the LIGO 1 and 
VIRGO detectors. It is therefore interesting to consider 
the distribution of the observed chirp masses as a potential 
observational statistic. One obvious advantage of such dis- 
tribution is that all the problems concerning the rate are 
suddenly vanishing. The uncertainties in the normaliza- 
tion that entered the calculation of the rate vanish when 
the distribution of observed chirp masses is considered. 
We shall now consider the sensitivity of the distribution 
of observed chirp masses to the cosmological model as- 
sumed, and to the sensitivity of a detector. Finally we 
will present the sensitivity of this statistic to the assumed 
model of stellar evolution. 

We will start by defining two reference models with 
which other models will be compared. These models cor- 
responds to the cosmological model described in section 
2, with the star formation rate given by the thick line in 
Figure 121 and the sensitivity of a detector given by equa- 
tion [D with Ai ^ 100 Mpc for model 1 and A2 = 1 Gpc 
for model 2 . In the following we will use a simple method 
of comparing distributions: the Kolmogorov Smirnov test. 
This test uses a parameters D defined as the maximum 
distance between two cumulative distributions. Two dis- 
tributions differing by D can be distinguished at a confi- 
dence level of approximately 10~^, when they are sampled 
at iV « 4/1)2 points. 

We first analyze the dependence of the expected distri- 
bution of observed chirp masses on detectors sensitivity. 
Variation in the sensitivity of a detector is modeled by 
the parameter Ai in equation ^ We present the results in 
Figure El for our two reference models. For the case of the 
reference model 1 (top panel of Figure , the maximum 
value of the parameter D reaches ~ 0.04 when the sen- 
sitivity is increased or decreased by a factor of two. For 
the more sensitive detector - reference model 2 in the bot- 
tom panel of Figure - the difference are larger and the 
parameter D reaches the value of ~ 0.15. This is due to 
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A [Mpc] 



Fig. 5. Variation of the shape of the distribution of ob- 
served chirp masses with the sensitivity of a detector. The 
top panel corresponds to the reference model 1, and the 
bottom panel to the reference model 2. On the vertical axis 
we plot the Kolmogorov Smirnov parameter D comparing 
a given distribution with the standard model described in 
the text. 

the fact that the redshift effects play a much stronger role 
for the more sensitive detectors. A smaller affects is con- 
nected with the fact that the detector starts seeing NSNS 
mergers from redshifts where the star formation rate was 
larger than locally. 

We present the dependence of the shape of the distri- 
bution of observed chirp masses on the assumed cosnio- 
logical model in Figures El and [3 The thick lines denote 
the case the reference model 1 with Ai = 100 Mpc, and 
we present the reference model 2 A2 — 1 Gpc with thin 
lines. Varying (while keeping the Universe flat, i.e. 
JIa = 1 — ^m) does not alter the shape of the observed 
distribution significantly: the parameter D never exceeds 
0.01 for model 2 and less than 0.003 for model 1. The as- 
sumed value of the Hubble constant has a stronger effect, 
however when Hq is varied by 20% the parameter D does 
not exceed 0.01 for the model 1 and 0.04 for the model 
2. We have also investigated the dependence of the distri- 
bution of observed chirp masses on the assumed shape of 
the star formation rate history. Apart from the standard 
model, with the star formation rate history given by the 
thick line in Figure 13 we investigated a model where the 
star formation rate ceases above z = 2, shown as a think 
line in Figure|3| and a model where the star formation rate 
is constant. The differences between the resulting distri- 
bution for the reference model and the model with star 
formation decreasing above z = 2 leads to D = 0.011, and 
for the flat star formation rate we obtained D = 0.021, 
while for model 2 these differences are: D = 0.010 and 
D = 0.036, respectively. 



Fig. 6. Change in the distribution of the observed chirp 
masses as a function of the assumed mass density in the 
cosmological model. The thick line corresponds to the ref- 
erence model 1, and the thin line to the reference model 
2. 



a 




40 60 

H [km s-i Mpc-i] 



Fig. 7. Change in the distribution of the observed chirp 
masses as a function of the assumed value of the Hubble 
constant. The thick line corresponds to the reference 
model 1, and the thin line to the reference model 2. 



We summarize the results of this section in Figures |S1 
and El Here we plot the values of the parameter D ob- 
tained from comparing our reference model 1 (Figure jSl , 
and 2 (Figure jSJ with distributions obtained when the 
stellar evolution is described by models listed in Table 1. 
We also present the maximal values of the parameter D 
obtained when we varied the detector sensitivity is var- 
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Fig. 8. The dependence of the parameter D comparing 
the distribution of observed chirp masses of the reference 
model 1 with other models of stellar evolution. For com- 
parison we also present the maximal values of D obtained 
when cosmological parameters are changed, for different 
star formation models, and when the sensitivity of the 
detector varies. 
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Fig. 9. Same as Figure|Hlbut for reference model 2. 



led, the cosmological model, and also for different assumed 
forms of star formation rate history. For most of the mod- 
els the value of D lies between 0.1 and unity. The excep- 
tion is model J for which D = 0.07. This is the model 
with different slope of the initial mass function. The pa- 
rameter D is small in this case since the distribution of 
observed chirp masses is dominated by heavy double black 
hole mergers. Decreasing the slope of the initial mass func- 



tion only enhances this effect, since it increases the ratio 
of black hole formation rate to neutron star formation 
rate. For a number of models: El, Gl, G2, and O the 
parameter D is larger than 0.5. These are the models 
which vary from the standard one by these parameters 
that affect the masses of compact object the most. The 
remaining models have the D parameter in the range be- 
tween 0.1 and 0.5. Distributions differing hy D = 0.2 can 
be distinguished when sampled with approximately one 
hundred points. Thus we c onfirm our previous conclusion 
I Bulik fc BelczvnskilEool that observation of about 100 
mergers shall yield significant constraints on the models 
of evolution of high mass stars. In Figure |H1 the expected 
distribution of chirp masses is mostly influenced by the 
underlying stellar evolution model, while in the case of 
reference model two shown in Figure variation of the 
sensitivity of the detector may influence the shape of the 
chirp mass distribution at the level similar to several stel- 
lar models. Thus we conclude that the distribution of ob- 
served chirp masses is a rather robust estimator of the 
underlying stellar evolution. 

5. Summary 

Using the Star Track binary population synthesis code we 
have investigated the properties of the population of dou- 
ble compact objects - the primary candidate sources of 
gravitational waves for high frequency interferometric de- 
tectors. We find that the distribution of lifetimes of various 
types of binaries is different: double neutron star binaries 
live typically for a few tens of million years, while the 
lifetimes of binaries containing black holes extend up to 
the Hubble time. Therefore, the double black hole binaries 
merging currently in our Milky way neighborhood origi- 
nate in systems formed even in the early star formation 
episodes! Hence, a proper calculation of the number and 
properties of the merging compact object binaries should 
include the full cosmological model of formation of these 
sources. We outline such calculation and calculate the ex- 
pected differential rate per unit observed chirp mass. We 
discuss the usefulness of the observed rate for determin- 
ing the properties of the underlying population. We es- 
timate the systematical uncertainty of the calculation of 
the rates due to uncertainty in the model. We find that 
the rate calculation carries a huge systematic error, com- 
parable to the spread in the rate due to variation of the 
stellar evolution model. We conclude that it is very un- 
likely that significant constraints can be obtained from 
consideration and modeling of the rate alone. However, 
we analyze the properties of another observable, the dis- 
tribution of observed chirp masses. We show that is very 
insensitive to the parameters that made the rate estimate 
uncertain. Moreover, we test the sensitivity of this statis- 
tic to the assumed parameters of the cosmological model: 
the mass density in the Universe, the value of the Hubble 
constant, and the model of star formation rate. All of these 
parameters hardly affect the shape of the distribution of 
observed chirp masses. We also verify that the shape of 
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this distribution does not vary significantly when a detec- 
tors sensitivity changes. This is especially important for 
detectors with non stationary noise. Non stationary noise 
would make comparison of the theoretical rate with ob- 
servations even more uncertain and difficult, yet it poses 
no problem for the analysis of the distribution of observed 
chirp masses. Analysis of the distributions of the observed 
chirp masses is therefore a very valuable tool for using 
the gravitational wave data to impose constraints on the 
stellar evolution models. 

The changes of the distribution of observed chirp 
masses with varying detectors sensitivity are mainly 
due to the fact that the observed quantity is the red- 
shifted chirp mass. Yet a redshift of a coalescing sources 
may be measured pro vided that we know its location 
l)Finn fc Chernof3ll99,'^ . A location can be estimated us- 
ing a network of gravitational wave detectors. A measure- 
ment of of the redshift would lead to the estimate of the 
chirp mass, and to removal of the bulk of the dependence 
of the observed chirp mass distribution on the fluctuations 
of detectors sensitivity. 

We must note that in general each measurement of a 
coalescence may also carry more information, like e.g. the 
individual masses of the coalescing object or their spins. 
Once these become available they must be included in 
the analysis. In this paper we have considered only the 
chirp mass measurement, as this is the most conservative 
approach. 
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